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SUMMARY 


The  functions  of  the  Keplerian  elliptic  motion  are 
expanded  with  respect  to  mean  anomaly  and  the  eccentricity 
by  applying  Poincare's  method  of  continuation  in  a  direct 
manner  no  the  equations  of  motion.  Two  algorithms  are 
proposed;  both  lead  to  programs  by  which  tne  classical 
expansions  are  constructed  symbolically  and  automatically 
on  a  computer  in  double  precision  arithmetic.  Even  to  as 
high  a  degree  as  30  in  the  eccentricity,  the  procedures  are 
remarkably  swift.  Special  care  is  taken  here  in  describing 
significant  error  controls  at  eacn  step  of  the  recurrences. 
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Astronomers  engaged  in  Celestial  Mechanics  devote  a  substantial 
portion  of  their  time  to  the  routine  manipulation  of  literal  expansions. 
It  has  long  been  recognized  that  digital  computers  are  capable  in 
principle  of  easing  this  burden  [3],  [4]. 

As  a  possible  start  toward  the  practical  implementation  of  that 
capability,  we  present  here  a  set  of  algorithms  which  have  been  written 
into  programs  for  an  IBM  7094.  The  package  deals  with  the  developments 
for  the  elliptic  motion  in  the  classical  Problem  of  Two  Bodies. 

The  mathematical  entities  to  be  processed  are  (truncated)  Fourier 
series  of  the  mean  anomaly  £,  their  coefficients  being  (truncated) 
power  series  of  the  eccentricity  e : 

(1)  F(£;e)  =  £  |" (  £  Cu  J  *  +  (  I  sk  ^sin  J ft"] . 

j>0  L'  k>0  /  Vk>p  '  -I 
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The  coefficients  C,  and  S,  .  are  rational  numbers,  and  it  would  be 

J  J 

in  keeping  with  the  ideal  spirit  of  symbolic  computing  to  store  and 

manipulate  them  as  such  in  the  machine.  Recently  [6],  this  formidable 

task  has  been  completed  with  success.  Nonetheless  the  newly  extended 

Cayley's  tables  are  not  likely  to  abate  the  discussion  whether  or  not,  in 

the  present  state  of  computing  machinery,  we  should  compromise  and 

represent  the  coefficients  C,  .  and  S,  .  as  floating  point  variables. 

K>J  K>J 


The  bolutions  of  the  Problem  of  Two  Bodies  are  known  in  a  closed 
form;  thus  their  analytical  expansions  in  d'Alembert  series  of  the  mean 
anomaly  and  the  eccentricity  serve  little  purpose  within  that  problem 
itself.  But  these  developments  are  the  indispensable  ingredients  of 
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elaborate  Perturbation  Theories.  Now  such  theories  operate 
spontaneously  in  the  field  of  real  numbers  and,  although  it  is 
often  of  the  utmost  importance  to  decide  how  close  the  multiple 
argument  of  a  trigonometric  term  is  to  a  commensurability  ratio, 
we  cannot  think  of  an  instance  where  it  was  essential  to  identify  the 
coefficient  itself  of  such  a  trigonometric  term  as  being  a  rational 
number.  Thus  we  find  some  ground  for  submitting  a  decisive  question: 
should  we  hold  for  more  than  a  mathematical  curiosity  the  fact  that 
the  coefficients  in  the  Fourier  expansions  of  a  Keplerian  motion  are 
genuinely  rational  numbers? 

Should  we  decide  at  the  beginning  of  a  Perturbation  Theory  to 
process  the  special  functions  of  the  Keplerian  motion  in  their  genuine 
arithmetic,  which  is  that  of  rational  numbers,  then  we  impose  on 
ourselves  the  formidable  task  of  operating  throughout  with  integer 
variables  of  multiple  length,  to  track  down  at  every  step  the 
nonrational  quantities  which  are  bound  to  occur  and  to  replace  them 
by  literal  symbols.  These  will  have  to  be  treated  by  the  program  as 
symbolic  variables;  in  this  way  the  result  will  involve  no  round  off 
error,  and  the  dependence  on  these  symbols  will  be  explicitly  displayed. 
But  thereby  the  number  of  auxiliary  variables  as  locumtenens  for 
irrational  numbers  will  have  proliferated  at  a  fast  pace.  Keeping  this 
growth  under  control  should  require  introducing  the  basic  concepts  of 
Galois  Theory  of  Fields,  especially  the  fundamental  distinction  between 
algebraic  and  transcendental  extensions  of  a  field.  Every  time  an 
irrational  symbol  is  introduced,  it  ought  to  be  decided  whether  it  is 
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transcendental  or  algebraic  over  the  field  of  rational  numbers,  and 
in  the  latter  case,  whether  or  not  it  belongs  to  the  algebraic 
extension  generated  by  the  algebraic  numbers  previously  introduced. 

Away  from  these  artificial  complexities,  our  aim  is  to  generate 
in  a  straightforward  manner  the  coefficients  in  double  precision 
arithmetic  for  the  expansions  of  the  elliptic  motion. 

At  the  suggestion  of  Moulton  [7],  [8],  development  for  the 
cosine  of  the  eccentric  anomaly  E  is  derived  immediate1  y  from  the 
differential  equation  for  the  radial  motion.  The  construction  consists 
mainly  in  applying  Poincare's  method  of  continuation:  a  transformation 

of  coordinates  exchanges  the  study  near  the  circular  solution  for  a 
study  near  an  equilibrium  position,  so  that  the  family  of  elliptic 
motions  having  a  preassigned  apsidal  line  and  a  given  mean  motion  turns 
out  to  be  the  family  of  periodic  orbits  issued  from  the  equilibrium 
point  representing  the  circular  trajectory  with  the  same  mean  motion. 

As  we  adopt  the  eccentricity  as  the  analytical  coordinate  along  the 
family  emanating  from  the  equilibrium,  the  method  of  analytical 
continuation  generates  the  expansion  of  cos  E  as  a  Fourier  series  of 
the  mean  anomaly  whose  coefficients  are  power  series  of  the  eccentricity. 
The  outstanding  merit  of  the  algorithm  is  that  it  generates  at  the  same 
time  the  powers  cos^E  as  series  of  the  same  type.  Consequently  only 
elementary  manipulations  of  i irnple  Poisson  series  are  necessary  to 
generate  therefrom  the  classical  expansions  in  the  Problem  of  Two  Bodies. 

As  an  alternative,  we  examined  another  algorithm  by  which  we  expand 
in  a  direct  manner  the  Cartesian  coordinates  x  ■  r  cos  f  and  y  ■  r  sin  f 


in  terms  of  the  mean  anomaly  and  the  eccentricity.  A  time  dependent 
transformation  of  coordinates  represents  the  circular  motion  as  an 
equilibrium  position  in  a  conservative  dynamical  system  with  two  degrees 
of  freedom.  Although  two  of  the  characteristic  exponents  at  the 
equilibrium  are  zero,  Poincare's  method  of  continuation  establishes 
that  the  equilibrium  generates  a  one-parameter  family  of  periodic 
solutions  which  actually  represents  the  elliptic  motions  having  a 
preassigned  line  of  apsides  and  the  same  mean  motion  as  the  generating 
circular  orbit.  Together  with  the  coordinates  x  and  y,  the 
analytical  continuation  produces  all  monomials  up  to  a  given 

degree  n  *  p  +  q. 

Both  direct  algorithms  proposed  here  call  upon  a  package  of 
subroutines  [2]  to  manipulate  series  of  the  type  (1)  Each  series 
is  denoted  by  a  FORTRAN  variable.  It  points  to  a  header  word  in  core 
which  contains  the  address  of  the  FORTRAN  variable  and  the  number  of 
terms  in  the  series.  The  miscellaneous  series  are  stored  dynamically 
in  sequential  order  in  a  preassigned  area  of  the  core  or  on  a  disk. 

As  it  was  first  proposed  by  Herget  and  Musen  [5],  a  series  is  considered 
as  a  mapping  F  of  the  product  set  N  x  N  *  Z/2  into  the  set  of  real 

numbers  (N  denotes  the  set  of  natural  integers,  Z  that  of  rational 
integers,  and  Z/2  the  set  of  congruences  modulo  2  in  Z  ,  thus  a 
set  reduced  to  the  two  elements  0  and  1.  Each  element  (j,k,I)  of 
N  x  N  x  Z/2  are  packed  in  one  word;  the  fact  that  F(j,k,I)  is  the 
image  of  the  element  (j,k,I)  by  the  mapping  F  is  represented 
physically  in  the  core  of  the  machine  by  placing  the  double  word 
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containing  the  real  number  F(j,k,I)  next  to  the  label  word  containing 
the  indices  (j,k,I).  The  coefficients  C.  ,  and  S.  ,  of  a  series 

j  »  K  J 

(1)  are  stored  below  its  header  word  dynamically  on  the  principle:  first 
come,  first  served.  I  -  0  refers  to  a  cosine  term,  I  ■  1  to  a  sine 
term;  the  multiple  j  of  the  argument  l  is  restricted  to  be  positive  for 
a  cosine  term  and  strictly  positive  for  a  sine  term  (by  changing  the 
sign  of  the  coefficient  S,  ,  ,  when  necessary). 

K  »  J 

To  automatize  completely  Poincare's  method  of  continuation,  we 
found  it  very  convenient  to  split  up  series  of  type  (1)  into  the  sums 


F(£;e)  -  J  F (i)ek 
k>0 


so  that  the  terms  appear  as  genuinely  finite  (and  not  truncated) 


j  i 

A  ’ 
\  ■ 


Fourier  sums  of  the  mran  anomaly  £.  Poincare's  continuation  is  then 

accomplished  not  by  successive  iterations  but  by  recurrence:  once  the 

coefficients  of  all  powers  of  e  have  been  determined  up  to  order  N, 

N+l 

that  of  e  is  computed  once  and  for  all.  At  this  stage,  the  key  of 

the  method  is  the  symbolic  inverse  of  a  differential  operator  to  be 
applied  to  a  trigonometric  sum.  But,  with  respect  to  the  natural  laws 
of  addition,  multiplication  and  scalar  multiplication,  the  series  of  the 
type  (1)  are  endowed  with  a  structure  of  commutative  graded  algebra  with 
a  unity  over  the  field  of  real  numbers;  moreover,  the  differential 
operators  which  come  up  at  various  stages  of  the  algorithms  are 
endomorphisms  if  not  of  the  whole  algebra,  at  least  of  some  of  its 
most  remarkable  subalgebras.  This  is  precisely  the  reason  why  Poincare's 
method  of  continuation  is  so  easily  implemented  automatically  on  a 
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computer. 
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1.  THE  DIRECT  EXPANSION  OF  cos  E 
In  the  Lagrangian  function 

(2)  £  -  Js(x2  +  y2)  +  ^  ,  r  -  |x2  +  y2|!*, 

which  describes  the  Problem  of  Two  Bodies  in  reference  to  a  fixed 
heliocentric  coordinate  system,  we  introduce  new  independent  and 
dependent  variables  by  means  of  the  definitions: 


(31)  x  -  aX, 

(32)  y  -  aY, 

(4)  n2a3  -  u. 

(5)  l  -  nt; 

we  also  define  a  normalized  radial  distance  R  such  that 
(61)  r  -  aR,  (62)  R  -  |X2  +  y2^, 

and  we  define  by  D  the  differentiation  with  respect  to  £.  Thus 
2  2 

on  dividing  by  a  /n  ,  the  transformed  Lagrangian  function  is 

(7)  £  -  (DX) 2  +  (DY)2]  +  j  . 

We  change  from  the  Cartesian  coordinates  (£,n)  to  the  polar 
coordinates  (R,0)  such  that 


(8) 


X  ■  R  cos  0, 


Y  ■  R  sin  0. 


_  _ _ — #.J 
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Since  6  is  ignorable  in  the  Lagrangian,  we  substitute  to  £  the 
Routhian  function  [10]: 


(9) 


where  C  is  the  constant  of  integration  introduced  by  the  integral 
of  angular  momentum 


(10) 


R  •  D0  -  C. 


Assume  that  the  conditions  for  an  elliptical  orbit  are  satisfied. 
Then  a  is  the  major  semi-axis  of  the  orbit,  n  the  mean  motion;  the 
constant  C  is  related  to  the  eccentricity  e  by  the  relation 


(ID 


C  -  vl-e  , 


and,  if  we  put  the  origin  of  time  at  an  instant  of  passage  at  perigee, 

£  is  the  mean  anomaly.  Denoting  by  E  the  eccentric  anomaly,  we  put 


(12) 


p  ■  cos  E, 


so  that  we  can  express  the  radius  vector  as  the  function 


(13) 


R  ■  1  -  pe. 


'*  v,-y«i rr  • 
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We  substitute  p  for  R  in  the  Routhian  function  (9),  and  we  expand 

2  2 
l/(l-pe)  and  l/(l-pe)  as  power  series  of  e.  On  dividing  by  e  , 

we  obtain  3t  in  the  form  of  the  power  series 

(14)  3-  *5<Dp)2  -  h  p2  +  \  J  (n  +  1)  (1  -  p2)pnen, 

n>l 

2 

where  we  have  omitted  the  constant  term  (1+e  )/2.  Thus  the  Lagrangian 
equation  becomes 

(15)  (D2  +l)p.h  (n  +  1) [n  -  (n+2)p2]pn'1en. 

n>  i 

Obviously,  when  e  *  0,  the  equilibrium  solution  p  =  0  represents 
the  Keplerian  circular  motion  with  radius  a  and  mean  motion  n.  Moreover 
the  variational  equation  around  that  equilibrium 

(16)  (D2  +  l)p  -  0 

admits  as  a  solution 

(17)  Pq  *  cos  i 

for  the  initial  conditions 

(18)  pQ(0)  -  1,  DP()(0)  =  0. 

Therefore,  according  to  Poincare's  method  of  continuation,  such  a  solution 
can  be  extended  analytically  into  the  power  series 


V 
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(19) 

whose  coefficients 

(a)  PnU  +  2tt  ) 

(b)  pn (0)  =  0, 


p(£;e)  -  £  Pn(*)en 

nto  n 

Pn(£)  for  n  >_  1  have  the  following  properties: 

pn^  ("condition  of  periodicity"); 

Dp  (0)  *  0  ("initial  conditions") . 
n 


On  substituting  (19)  in  the  original  differential  equation,  we 
obtain  from  equating  like  powers  of  e  in  both  members  a  recurrent 
system  of  differential  equations  which  determine  the  coefficients  of 
(19).  Before  we  write  down  this  system,  let  us  introduce  a  convenient 
notation.  We  put 


(20)  Pn(fc;e)  *  I  p^n)a)eP  (n  >  0); 

p>0  P 

thus  p^  represents  the  Fourier  series  which  is  the  coefficient  of 
eP  in  the  expansion  of  (cos  E)n.  With  this  convention,  we  can  write 
the  coefficient  of  en  in  the  right-hand  member  of  (15)  as  the  function 


(21)  UnU) 


1 

2 


2 

l<k<n 


(k  +  1)  kp 


(k-1) 

n-k 


-  (k  +  2)p 


(k+1)' 

n-k 


(n  >  1) , 


and  we  find  that  p 

n 


is  a  solution  of  the  differential  equation 


(22) 


(D2  +  l)pn 


U  . 
n 


The  following  facts  can  be  established  by  recurrence  on  the  index  n 
from  examining  the  right-hand  members  and  taking  into  account  the 

initial  conditions: 
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(a)  For  any  n  ^  0,  is  a  finite  trigonometric  sum;  more 

precisely,  it  does  not  contain  multiples  j£  of  the  mean  anomaly 

beyond  j  -  n  +  1. 

(b)  For  any  n  >  0,  p  is  an  even  function  of  £ ;  thus  it 

—  n 

contains  only  cosine  terms. 

(c)  For  any  even  n  >  0,  p  (£  +  tt)  ■  -p  (£) ,  whereas  for  any 

—  n  n 

odd  n  >  1,  p  (£  +  tt)  =  p  (£). 

—  n  n 

Property  (a)  guarantees  that  the  coefficients  p^  can  be  computed 
from  the  recurrence  in  a  closed  form.  In  other  words,  they  do  not  result 
from  successive  approximations,  but  once  their  respective  differential 
equation  (22)  has  been  solved,  they  are  known  as  accurately  as  it  can  be, 
save  for  errors  resulting  from  cumulative  rounding  off  in  the  arithmetic 
manipulations  of  the  coefficients. 

Property  (c)  implies  that,  at  an  even  order  N  =  2m,  the  right-hand 
member  contains  only  cosine  terms  of  odd  multiples  of  the  mean 

anomaly.  Thus  it  looks  as  if  it  would  contain  among  others  a  term  in 
cos  £  .  As  this  would  prevent  any  solution  of  (22)  from  being  periodic, 
we  have  to  interpret  Poincard's  theorem  as  meaning  that,  at  every  even 
order,  the  various  contributions  to  arrange  themselves  so  as  to 

cancel  exactly  the  term  in  cos  £.  In  the  course  of  the  recurrence, 
round  off  errors  will  prevent  an  exact  cancellation.  Thus  the  coefficient 
of  cos  £  in  will  provide  a  control  on  the  errors  at  every  even 

order  of  the  recurrence. 

We  list  in  Table  I  the  size  of  the  coefficient  of  cos  £  as  it 
appeared  in  the  right-hand  members  of  Uj  before  these  spurious  terms 
were  rejected  by  the  program. 


At  the  odd  orders,  the  recurrence  brings  two  controls  on  the 
error. 


1°)  Let  us  recall  that  the  average  value  of  r/a  over  the  mean 

2 

anomaly  is  equal  to  1  +  e  /2.  Consequently,  the  coefficients 
P2n+1  ^  should  be  exempt  from  constant  terms,  which  implies 

that  the  right-hand  members  l^n+l  themselves  should  be  exempt 
from  such  terms.  In  point  of  fact,  in  the  course  of  building 
^2n+l*  they  cancel  out  exactly  when  the  various  pn  and  their 
powers  are  produced  in  integer  arithmetic.  But,  in  our  algorithm, 
as  we  proceed  in  double  precision  arithmetic,  they  might  still 
appear,  but  as  very  small  numbers.  Hence  they  serve  as  tests  on 
the  adverse  effects  of  round  off  errors. 

2°)  In  order  to  integrate  the  differential  equations  (22),  we 
proceed  in  two  steps.  First  we  geneiate  the  particular  solution 


4 

i 


J 
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(23) 


d2+i 


n 


(n  >  1) 


2 

by  application  of  the  inverse  differential  operator  1/D  +1 
according  to  the  rules 


(20  cos  p £  -  — 5 —  cos  pi,  (p  +  1), 

D  +1  -p  +1 

(24?)  sin  p £  =  — ~ —  sin  pi,  (p  4  1). 

D  +1  -p  +1 

As  we  just  observed,  this  operation  is  always  possible:  neither  at 

odd  nor  at  even  orders,  does  U  contain  critical  terms  in  cos  £. 

n 

Then  to  the  particular  solution  p*,  we  should  add  the  general 
solution  of  the  homogeneous  differential  equation,  so  as  to  obtain 


p  =  C  cos  £  +  S  sin  £  +  p*, 
n  n  n  n 

where  the  arbitrary  coefficients  and  ought  to  be  determined 

by  the  initial  conditions. 

Of  course,  since  p  is  a  sum  of  cosine  terms  only,  we  find  at  any 

n 

order  that  S  =0. 

n 

At  the  even  orders,  we  find  that 


C2n=-pSn(°>- 


> 
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But  at  the  odd  orders,  we  know  that  P2n+i  and  U2n+1*  hence  P$n+1  ’ 
contain  only  cosines  of  even  multiples  of  l.  Therefore,  not  only  should 
we  make  ^2n+^  “  0,  but  we  should  also  check  that  the  particular  solution 
C*  meets  the  initial  condition 


In  Table  II,  we  enter,  for  each  odd  order,  the  size  of  the  spurious 


independent  term  which  came  into  ^n+l  anc*  t*ie  values 


1n+l(0>- 


Table  II 
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2.  DERIVED  EXPANSIONS 

Other  special  functions  of  the  Keplerian  elliptic  motion  may 
be  expanded  by  elementary  algebraic  operations  from  the  developments 
of  cospE  (p  >_  1)  that  we  have  just  obtained. 

Thus,  from  the  series 


cosPE  =  pP  =  ][  PnP^n  (PiD. 

n>0  n 

we  find  immediately  the  expansion  of 


(25) 


R  =  ~  =  1  -  e  cos  E  ~  1  -  V  p(1J  e", 
a  Kn-1 

n>l 


(26) 


1  a  _  1 


R  r  1-e  cos  E 


v  n  n  v 

I  p  e  =  I 

n>0  n>0 


(P) 


l  p'P 

rv.  n-p 


p>0 


(27)  x  =  *  =  e  -  cos  E  =  -1  +  (1  -  )e  -  ^  p| 

n>2  n 


(l)en 


(28) 


E  =  [  - 

J o  r 


(s) 


ds . 


We  observe  that  the  average  of  a/r  over  the  mean  anomaly  is  equal  to  1, 
so  that  E  -  Z  turns  out  to  be  a  d'Alembert  series  in  l  and  e  whose 
average  over  the  mean  anomaly  is  zero.  Therefore,  by  elementary 
manipulations  of  series  of  the  type  (1),  we  are  able  to  expand  in  the  mean 
anomaly  the  function 


sin  E  =  ~  (E  -  £) , 


(29) 
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and  this  development  in  turn  yields  the  expansion  of 

(30)  Y  *  1  -  '/l-e2  sin  E. 

a 

The  "equation  of  the  center"  is  given  by  the  quadrature 

f  =  4-e2  /  —  ds. 

J0  K  (s) 

To  experiment  more  extensively  with  the  simple  programming  techniques 
proposed  here,  we  undertook  to  expand  all  powers  (r/a)P,  2  p  15,  in 
series  of  the  mean  anomaly  and  of  the  eccentricity.  The  computation  is 
easy,  since  it  reduces  to  repeated  additions  and  multiplications  of 
homogeneous  components.  But  as  the  exponent  p  increases,  the  round  off 
errors  get  bigger  and  bigger.  We  have  two  ways  of  watching  the  increasingly 
destructive  effects  of  the  round  off  errors. 

First  we  know  that,  for  any  integer  p  ^  0,  the  term  of  either 

(r/a)2p  1  or  (r/a)2p  which  is  independent  of  l  is  a  polynomial  in  e2 

whose  degree  is  2p.  Now  the  program  ignores  this  property,  and 

consequently  this  polynomial  appears  in  the  output  as  a  (truncated)  series 
2 

in  e  .  We  can  estimate  the  contaminating  influence  of  the  round  off  errors 
by  inspecting  these  spurious  coefficients.  They  are  entered  in  Tables  III 
and  IV. 

Another  check  on  the  calculations  of  the  successive  powers  (r/a)P 
is  suggested  by  Brown  and  Shook  [1].  For  any  rational  integer 
p  (>  0  or  <  0) ,  we  find  that 


us,** 
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Table  IV.  Error  control  on  the  Independent  terms  of  (— )  (continued) 
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DRP  -  pRP_l  •  DR, 

D2RP  -  p(p  -  1)RP'2  *  (DR)2  +  pRP_1  •  D2R, 
so  Chat,  In  view  of  the  energy  integral 

2 

(31)  (DR)2  +  -  ~  =  -1 


belonging  to  the  Routhian  function  (9)  and  of  the  Lagrangian  equation 


(32) 


derived  from  the  same  function  (9),  we  obtain  that 


(33)  -jj  D2Rp  +  (p  -  1)RP  2  -  (2p  -  3)RP“3  +  (p  -  2)  (1  -  e2)RP~4  =  0. 


Thus,  for  instance,  by  making  p  equal  to  2,  3,  4  successively. 


■j  D2R3  -  -2R  +  3  -  (1  -  e2)  £  , 
£  D2RA  =  -3R2  +  5R  -  2(1  -  e2). 


The  left-hand  members  of  (33)  are  not  as  good  an  error  control  as  the 
residuals  in  the  power  series  of  e  which  is  the  independent  term  of 
(r/a)P.  Indeed  they  imply  that  we  take  the  second  derivative  of  a 
Fourier  series,  and  that  is  not,  as  we  know,  the  most  favorable  operation 
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we  should  perform  on  a  trigonometric  sum.  We  do  not  reproduce  here 
the  series  which  we  have  computed  from  (33)  for  p  *  5  to  20;  their 
coefficients  are  somewhat  similar  to  the  residuals  listed  in  Tables  III 
and  IV. 

On  the  whole,  a  detailed  examination  of  the  round  off  errors 
produced  in  any  coefficient  of  the  powers  (r/a)11  and  (a/r)^ 

(1  p  £  20)  has  convinced  us  that  the  elementary  techniques  advocated 
here  guarantee  more  than  the  accuracy  which  is  required  in  Perturbation 
Theories  from  the  expansions  of  the  elliptic  motion. 

Other  techniques,  of  course,  are  possible.  But,  as  we  shall  show 
it  on  another  algorithm  in  the  next  section,  it  is  of  the  utmost 
importance  that  Jhe  residuals  should  be  examined  carefully  before  a 
proposed  technique  is  adopted,  and  its  re*,  i  Its  entered  as  data  into  a 
major  Perturbation  Theory. 

3.  THE  DIRECT  EXPANSION  OF  (r/a)  cos  f  AND  (r/a)  sin  f 

Let  us  come  back  to  the  Lagrangian  function  (7)  expressed  in 
normalized  heliocentric  Cartesian  coordinates  (X,Y).  To  these  variables 
we  substitute  the  variables  (£,n)  defined  by  the  non-conservative 
relations 


(34) 


X  - 

cos 

l 

+ 

e 

cos 

l  - 

n 

sin 

l. 

Y  - 

sin 

l 

+ 

sin 

l  + 

n 

cos 

l* 

so  that  the  Lagrangian  function  becomes 
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(35)  £  -  *5[(D02  +  (Dn)2]  +  a  •  Dn  -  n  •  DO  +  ^R2  +  j 
with 

(36)  R  =  |l  +  2£  +  C2  +  n2  j^. 

It  is  obvious  that  £  =  n  =  0  is  an  equilibrium  solution  for  the  system; 
in  fact,  it  corresponds  to  the  circular  motion.  In  order  to  study  the 
Keplerian  motions  around  this  equilibrium,  we  expand  R  in  Taylor  series 
in  the  powers  of  £  and  n- 

Such  an  expansion  is  carried  on  automatically  on  the  computer  by  means 
of  programs  which  manipulate  symbolically  double  power  series. 

We  look  at  a  series 


A(£  ,n) 


I  I 

j>0  k>0 


_j  k 

V  n 


as  a  sum  of  homogeneous  components 


A  -  £  A^,n), 

n>0 


where,  for  any  n  ^  0, 


A 

n 


I 

0<j  <n 


n-j  ,  j 


Cn~V 


is  a  homogeneous  polynomial  of  degree  n  in  the  variables  £  and  n. 
2 

If  we  take  R  for  the  series  A,  we  would  have 
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Ax  “  2C, 

2  2 
A2  -  C  +  n  . 


Since  Aq  >  0,  the  series  A  has  a  formal  square  root,  which  we 

shall  denote  by  B.  In  our  case,  B  is  the  Taylor  expansion  of  the 

■A  _ 

normalized  distance  R  in  the  powers  of  £  and  n*  The  homogeneous 
components  B^  of  B  are  determined  by  recurrence  from  the  fundamental 
relation 


B 


2 


=  A 


according  to  the  formulae: 


B 


0 


A1/2B0’ 


B 

n 


2 

1< j  <n-l 


B.B  4 
J  n-j 


/2B, 


(n  12). 


Since  Bq  =>  1,  the  series  B  has  a  formal  inverse,  which  we  shall 
denote  by  C.  In  our  case,  C  is  the  Taylor  expansion  of  1/R  in  the 
powers  of  £  and  n.  The  homogeneous  components  of  C  are  computed 

recursively  from  the  definition 


BC  -  1 


according  to  the  formulae 


C 


0 


-B1C0/B0, 


C 

n 


2 

0<j  <n-l 


B  i c  j 

n-j  j 


/B, 


(n  >  2). 
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Through  these  manipulations  of  polynomials,  we  come  to  the 
development 

i  -  1  -  t  +  C2  -  hr?  +  I  (  I  /-1  A 

n>3  yO<j<n  nJ,J  j 

Accordingly,  around  the  equilibrium,  the  Lagrangian  function  £  is 
decomposed  into  the  series 


£  =  y  £  , 

n>2  n 


£2  =  ls[(D02  +  (Dr,)2]  +  U  •  Dn  -  n  •  DO  +  |  £2, 


£n  =  y  w  .  .£n"Jnj  (n  >  3). 
n  r\  ±  n— 1 ,  l  — 

0< j  <n  J  J 


We  address  ourselves  first  to  the  Lagrangian  equations  derived 
from  £ namely 


(37) 


(D  -  3)£  -  2Dn  =  0, 
D(Dn  +  20  =  0. 


They  constitute  a  homogeneous  linear  system  with  constant  coefficients. 
Using  Putzer's  algorithm  [9],  we  find  its  resolvent  R(£).  It  is  a 
4x4  matrix 


RU) 


(ri,jU))l<i,j<4 


I 


whose  coefficients  are  as  follows 


4-3  cos  £, 


-6 (l  -  sin  £)  , 


rl,2  =  °* 

r2 ,2  "  lf 

r1  3  =  sin  £, 

r2,3  *  ~2^  "  cos  * 

rlt4  =  2(1  -  cos  £)  , 

r2,4  =  +  4  sin  *. 

r3  1  *  3  sin  £, 

r4  “  -6(1  -  cos  £) , 

r3,2  =  °’ 

r4,2  =  ° 

r3,3  =  cos  *• 

r4>3  =  -2  sin  £, 

r3>4  =  2  sin  £, 

r4  4  =  ~3  +  4  cos  £. 

The  general  solution  of  (37) 


C  rl,lC0  +  rl,2n0  +  rl,3^0  +  rl,4n0’ 


n  =  r2,lC0  +  r2,2n0  r2,3C0  +  r2,4n0’ 


can  thus  be  written  in  the  form 

£  ■  2A  -  B  cos  £  +  C  sin  £, 

(38) 

n  =  D  -  3A£  +  2C  cos  £  +  2B  sin  £, 

where  the  arbitrary  constants  A,  B,  C,  D  will  have  to  be  determined  by 
the  initial  conditions. 


-*•.'**& * 
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In  view  of  (38),  the  first  order  variational  equations  around  the 
equilibrium  admit  a  periodic  solution  if  and  only  if  A  =  0.  Moreover, 
this  solution  (£^(4) ,n^(4))  satisfy  the  symmetry  conditions 

q(-0  =  SjU),  n1(-£)  =  -n1(£) 

if  and  only  if  C  =  D  =  0.  Hence  B  remains  as  a  parameter,  and  from 
what  we  know  of  the  Keplerian  motions,  we  can  identify  B  to  the 
eccentricity  e. 

To  summarize,  a  first  order  analysis  of  the  motion  around  the 
equilibrium  yields  a  one  parameter  family  of  symmetric  periodic  orbits 

=  -e  cos  4, 
n1  *  2e  sin  4. 

The  problem  is  to  continue  analytically  this  family  into  a  solution 
of  the  complete  differential  equations,  which  is  of  the  form 

(39)  £(4;e)  =  ^  £,  (4)ek,  n(4;e)  =  ]>  n.(4)ek 

k>l  R  k>l  K 

where  the  coefficients  and  are  periodic  in  4  with  period  2tt. 

In  view  of  the  symmetry  R(£,-n)  =  R(£,n)»  the  coefficients  w.  , 

J 

in  the  homogeneous  components  ^^(n  i  3)  vanish  for  all  odd  indices  k. 
Hence  we  can  impose  on  the  coefficients  £^(<0  and  h^(4)  the  symmetry 


conditions 
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(40)  Sk(-£)  =  \(-)0  “  -nkU)  (k  -  1^’ 

It  implies  that,  for  any  k  ^  1,  £k  is  a  sum  of  cosine  terms  only, 
and  nk  a  sum  of  sine  terms  only. 

We  shall  also  impose  the  perigee  condition 


(41) 


5k(0)  =0  (k  >  2). 


Using  Poincare's  method  of  the  small  parameter,  it  can  be  established 
that  there  exists  one,  and  only  one,  family  of  periodic  solutions  in  the 
form  (39)  which  satisfies  the  symmetry  conditions  (40)  and  the  initial 
conditions  (41). 

In  order  to  outline  the  recursive  steps  that  we  have  to  take  in 
determining  the  coefficients  and  n j ,  we  shall  introduce  the 

following  notations: 

1  k 

1°  For  any  j  >_  1,  and  k  >_  1 ,  we  write  the  monomial  £Jn  as 
the  series 


S 


V 


2  X;|’k(t>en. 

n>j+k 


2°  We  put 


U-  s  n 


n>3 


as  ’ 


d£ 

v  «  s  — a  • 


n>3 


these  partial  derivatives,  when  expressed  in  terms  of  e  and  l , 
are  series  of  the  form 
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I 

♦ 


U  -  I  Un(Oen,  V  =  X  VnU)en, 


n>2 


n>2 


wherein 


U  -  X  X  (k  +  Du...  ,Xk,m_k(£) 

n  „  *4.  k+1  ,m-k  n 

2<m<n  O'kcm 


V  -  y  T  <k  +  Du  ..,x""k,k(«) 
n  2<m<n  0<k<m  "-k'k+1  n 


With  these  notations,  it  is  obvious  that,  once  all  coefficients 

X^’  (£)  are  known  up  to  order  n  -  1  for  j  >_  1 ,  k  >_  1  and 

j  +  k  £  n  -  1,  then  we  can  generate  by  multiplication  of  finite 

i  k 

trigonometric  sums  the  coefficients  ’  (£)  of  order  n  for  j  >_  1, 
k  ^  1  and  j  +  k  =  n,  and  consequently  the  expressions  and  V  . 

In  this  way  we  have  built  the  differential  equations 


(42) 


(D  -  3K  -  2Dn  =  U  , 
n  n  n 


D(Dn  +  20  =  V 
n  n  n 


which  define  the  terms  of  order  n  in  the  scries  (39). 

Before  we  integrate  this  system,  let  us  state  two  properties  which 

are  satisfied  by  the  coefficients  i  and  n  • 

n  n 

(a)  For  any  n  >_  0, 


(43)  {2n+1U  +  .)  -  -t2„+1(t>,  n2n+1U  +  ">  -  -"2n+l(l) 


V 


In  other  words,  the  coefficients  at  odd  orders  contain  only  odd 
multiples  of  the  argument  l. 

(b)  For  any  n  >  1, 


(44) 


^2n(£ 


+  i0  =  C2n(£)’ 


n2n(£ 


+  it) 


n2n(0. 


At  even  orders,  the  coefficients  contain  only  even  multiples  of  the 
argument  l. 

These  properties  can  be  established  by  recurrence.  They  are  essential 
to  the  success  of  the  recurrence  scheme  by  which  the  coefficients  £  and 
nn  are  computed  and,  as  we  shall  point  out,  they  provide  valuable  checks 
on  the  accuracy  of  the  calculation. 

We  determine  a  particular  solution  of  the  system  (42)  by  first 
performing  the  quadrature  which  is  the  second  equation  and  by  substituting 
the  result  into  the  first  equation.  Thus,  if 


W 

n 


—  V 

D  V 


<J>  *  U  -  2W  (n  >  2)  , 

n  n  n  — 


we  find  that  the  system  is  reduced  to  the  equation 

(43)  (D2  +  1)E  =  ♦. 

n  n 


At  this  point,  we  ought  to  remark  that,  at  even  orders,  is  a  sum  of 

cosines  of  even  multiples  of  l,  but,  at  odd  orders,  it  is  a  sum  of  cosines 


of  odd  multiples  of  £.  That,  in  ;his  latter  case,  it  does  not  contain 
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i 

i 


l 


terms  in  cos  £  is  a  consequence  of  Poincare's  theorem.  In  point  of 
fact,  should  we  carry  the  calculation  of  tf^n+l  *n  *nte8er  arithmetic, — 
and  we  did  it  by  hand  up  to  order  5  — ,  we  would  find  that  this  critical 
term  vanishes  by  cancellation.  But,  as  we  compute  in  double  precision 
arithmetic,  it  will  appear  with  a  small  coefficient  about  the  size  of  the 
cumulative  errors  caused  by  rounding  off.  As  it  constitutes  a  valuable 
control  on  the  errors,  we  have  reproduced  it  in  Table  V. 

We  now  apply  the  rules  (24^)  and  to  construct  the  particular 

solution 


(44) 


d2+i 


4>  . 
n 


Observe  that  E*  contains  a  term  independent  of  £  if  and  only  if  n 
is  even.  Thus,  in  the  quadrature 


(45) 


"S  '  D(W„ 


2CJ) 


resulting  from  the  second  of  the  equations  (42),  there  occurs  a  secular 
term,  say  Kn£,  if  and  only  if  n  is  even. 

Because  the  equations  (37)  constitute  the  homogeneous  system  belonging 
to  (42),  the  general  solution  of  (42)  is  of  the  form 


E  *  2A  -  B  cos  £  +  C  sin  £  +  E** 
n  n  n  n  n 

n  =  D  -  3A  £  +  2C  cos  £  +  2B  sin  £  +  n*. 
n  n  n  n  n  n 

In  view  of  the  symmetry  conditions  (40),  we  have  that 


V 


i 


-29- 


C  -  D  -  0. 
n  n 

The  other  two  coefficients  A  and  B  are  to  be  determined  in  different 

n  n 

wavs  depending  on  whether  the  order  n  is  even  or  odd. 

(a)  If  n  is  even,  the  parity  rule  (44)  implies  that 

B  =  0; 
n 

and,  for  nn  to  be  periodic  in  £,  we  must  have  that 

A  =  -K  /3. 
n  n 

In  which  case,  the  perigee  condition  (41)  provides  a  check  on  the 
accuracy,  namely  the  relation 


1  - 


(46) 


2A  -  C*(0)  =  0. 
n  n 


The  residuals  on  this  relation  are  entered  in  Table  V. 

(b)  If  n  is  odd,  because  n*  is  exempt  from  secular  terms,  we 
must  have 


for  nn  to  be  periodic.  Then  the  perigee  condition  (41)  serves  to 
compute  that 


B 

n 


c;<o). 


l 
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Table  V 


Error  controls  in  the  expansions  of  £  and  n 


Order 

Error 

Order 

Error 

2 

0 

3 

0 

4 

0.6 

X 

io-15 

5 

0 

6 

0.4 

X 

H-* 

o 

1 

7 

0.1 

X 

io-13 

8 

0.4 

X 

io'13 

9 

0.7 

X 

10-13 

10 

-0.1 

X 

icf12 

11 

0.4 

X 

io-12 

12 

-0.1 

X 

io"11 * * V * * * * X 

13 

0.1 

X 

o 

o 

r—H 

14 

0.1 

X 

io-10 

15 

0.1 

X 

10  9 

16 

-0.1 

X 

io-9 

17 

-0.2 

X 

00 

1 

o 

H 

Compared  with  the  corresponding  ones  in  Tables  1  and  II,  the  entries 

of  Table  V  indicate  that  the  algorithm  we  just  described  is  less  favorable 

than  the  first  one  in  accuracy.  The  reason  is  that  the  series  and 

V  are  obtained  by  combining  terms  of  the  same  size — a  defect  already 

n 

mentioned  by  Brown  and  Shook  [1], 

Before  we  conclude  the  analysis  of  this  algorithm,  let  us  consider 

one  more  check.  Once  the  expansions  of  £  and  n  are  known,  those  for 

X  =  x/a  and  Y  =  y/a  can  be  calculated  from  the  transformation  formulae 
(34).  Now  we  know  that  the  average  value  of  X  over  the  mean  anomaly  is 
equal  to  -3/2  e,  which  means  that  the  independent  term  in  the  expansion 
of  X  reduces  to  that  term.  Actually,  by  composing  the  developments  of 
C  and  n  according  to  the  first  of  the  transformation  formulae  (34), 
we  find  for  the  independent  term  a  truncated  series  in  e.  In  the 
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second  column  of  Table  VI,  we  enter  those  coefficients  as  ihey  have  been 

computed  by  the  program;  they  reflect  the  increasingly  damaging  propagation 

of  the  round  off  errors  throughout  the  recurrence.  Indeed,  that  the 
average  value  of  X  does  not  contain  the  eccentricity  beyond  the  first 

power  implies  that,  for  each  n  >_  1,  the  coefficient  ^n+l  ^  of 

2n+i 

e  cos  £  in  the  expansion  of  £  should  be  equal  to  the  coefficient 
^2n+l  1  e^n+^sin  £  in  the  expansion  of  n.  For  the  sake  of 

comparison,  we  enter  the  difference  l^2n+l  1  ~  ^2n+l  1^  secon(* 

column  of  Table  VI. 


Table  VI 


Error  controls  in  the  expansion  of  x/a 


Order 


x/a 


lX2n+l.l~Y2n-H.l 


3  0  0 


5 

-0.4 

X 

10-15 

0.75 

X 

io-15 

7 

-0.1 

X 

io"13 

0.20 

X 

io-13 

9 

-0.4 

X 

10-13 

0.74 

X 

io-13 

11 

0.8 

X 

io-12 

0.16 

X 

io-11 

13 

-0.3 

X 

io-11 

0.56 

X 

io~n 

15 

-0.3 

X 

o 

1 

t~* 

o 

0.55 

X 

io-10 

17 

0.8 

X 

Os 

1 

o 

»— 4 

0.16 

X 

H-* 

o 

1 

00 

*<3***,.  as** . 
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4.  CONCLUSIONS 

The  classical  expansions  of  the  elliptic  motion  in  terms  of  the 
mean  anomaly  can  be  obtained  in  a  straightforward  manner  by  applying 
Poincare's  method  of  continuation  to  finding  periodic  solutions  to  the 
differential  motions  in  the  neighborhood  of  the  circular  motions. 

The  programs  that  we  propose  here  are  fast.  The  first  algorithm 
builds  the  expansions  of  all  powers  (cos  E)*3,  1  <  p  <_  30,  up  to  the 

degree  30  in  the  eccentricity  in  less  than  6  minutes  on  an  IBM  7094. 
Because  it  needs  a  disk  file  to  store  the  series  representing  the 
monomials  (0  P  1.  J7,  0  <_  q  <  17),  the  second  algorithm  takes 

18  minutes  of  an  IBM  7094  to  complete  the  recurrence  at  order  17. 

As  long  as  the  order  of  the  expansions  is  not  pushed  unrealistically 
too  high  up,  these  elementary  programs  can  easily  become  parts  of  major 
undertakings  as  required  by  the  classical  Perturbation  Theories. 

A  special  effort  was  made  to  provide  enough  error  controls,  so  that 
the  reader  may  judge  the  quality  of  accuracy  he  may  expect  from  these 
algorithms  at  various  orders.  Thus  we  concede  that,  although  it  is  more 
extensive  in  the  intermediary  results  that  it  brings  out,  the  second 
algorithm  proposed  here  yields  sensibly  less  accuracy  at  somewhat  higher 


orders  than  the  first  one. 
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ANNEXES 

To  the  reader  who  would  like  to  set  up  subroutines  needed  for 
manipulating  Poisson  series  of  the  type  (1,1)  and  to  assemble 
therefrom  a  program  to  apply  Poincare's  method  of  continuation  to  the 
Problem  of  Two  Bodies,  it  may  be  of  some  use  to  check  the  results 
generated  by  the  computer  against  the  correct  answers  at  the  few  first 
orders.  Thus  we  reproduce  here  the  calculations  made  by  hand  against 
which  we  checked  our  procedures.  The  notations  in  Annexes  I  and  II 
adhere  strictly  to  those  introduced  respectively  in  Sections  1  and  3. 
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Order  0 


Order  1 


Order  2 


Order  3 


ANNEX  I 

Expansion  of  cos  E 


Pq  =  cos  2. 

b'l  =  -  |  -  |  cos  2  2, 

P1  =  "  i  +  2  COS  2*-' 

ct  =  0, 

1,1 

P^  =  -  ~  +  cos  2  2. 


U2  =  -3  cos  32, 

3 

P*  =  g  cos  32, 

C  =  -  1 
2  8  ’ 

3  3 

i>2  =  -  g  cos  2  +  -  cos  3C. 

=  cos  22-5  cos  42, 

P*  =  ~  3^  cos  ^  ^  +  y  cos  42, 

C3  =  °’ 

P^  -  -  ~  cos  2  2  +  -j  cos  42. 


-39- 


Order  4 


Order  5 


»? 


45  ,,  125  ,, 

—  cos  3 2  -  ~rr  cos  5  2,, 
lb  lb 

45  x  125  .. 

"  128  C0S  3*  +  384  C0S  5£| 


C,  = 


P/.  = 


„  45  .  125 

cos  2  -  Y^g  cos  32  +  cos  51. 


U5  = 
P*  - 


3  189 

-  7T  cos  22  +  6  cos  4 2  -  ~r—  cos  61, 

lb  lb 

1  „  2  ,  .  ,  189  ,  . 

—  cos  21  -  —  cos  4 2  +  777T  cos  62, 

lb  d  obO 


0, 


1  on  2  /a  *  189 

pc  =  —  cos  21  -  —  cos  42  +  ~  cos  62. 

D  16  j  DbU 


These  coefficients  were  checked  by  substituting  the  solution  in  the 
energy  integral  belonging  to  the  Routhian  function  JP  given  by  formula 


(9)  in  the  text. 


-40- 


Order  1 


Order  2 


Order  3 


ANNEX  II 

Expansion  of  £  and  n 


“  -e  cos  4, 

=  2e  sin  4. 

2 

Coefficient  of  e 

2  2  2 

in  C  '  he  +  he  cos  24, 

2 

in  £n :  -e  sin  24, 

2  2  2 

in  n  J  2e  -  2e  cos  24; 

3 

U2  -  J  cos  21 • 

3 

V2  -  J  cos  2l* 

Z*  =  2  +  2  COS  2l’ 

-  -34+  h  sin  24, 
a2  =  —1,  B2  *  0, 

C2  *  ~  "I  +  2  COS  2<l’ 

n2  *  h  sin  24. 

Coefficient  of  e^ 

2 

in  C  :  *s  cos  4  -  h  cos  34, 

13  3 

in  - g"  sin  4  +  -g-  sin  34, 

2 

in  n  :  h  cos  4  -  h  cos  34, 
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Order  4 


in 

C3: 

in 

£2n 

in 

Cn2 

in 

3. 

3 

4 


1 

4 


U„  = 


V.  = 


in  £  n:  *5  sin  £  +  sin  32., 


6  sin  2  -  2  sin  32; 

25  TO 

~r  cos  3?., 


■7  cos  2 
4 


9  39 

g-  sin  2  — g-  sin  32 , 


=  Q  C0S  3*. 

9  7 

=  -  g-  sin  2  +  —  sin  32, 


A3  -  °,  »3  -  g  , 

3  3 

=  -  g  cos  2  +  g-  cos  32, 

3  7 

r) 3  =  -  —  sin  2  +  ^  sin  32. 


Coefficient  of  e 


in 

C2: 

3  1 

4  -  2  C°S 

22  - 

1  /  0 
cos  42, 

in 

Cn: 

-  1  sin  22 

0 

sin  42, 

in 

2 

n  : 

23  .  4  . 

~  32*  +  *3  cos 

-  gg  cos  42 

in 

C3: 

3  .  3 

-  g  +  g  cos 

*2, 

in 

£2n: 

9  ■  TO 

:  —  sin  22 

7 

16 

sin  42, 

in 

f,n2: 

:  -  7  +  2  cos  2£ 
4 

-  —  cos  42, 
4 

in 

3 

n  : 

|  sin  22  - 

*7  sin  4£ , 

4 

in 

54-. 

3  ^  1 
•g  +  —  cos 

22  + 

1  /  0 

-  cos  42, 
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Order 


3 

in  £  n:  sin  2i  -  h.  sin  4 1, 

2  2 

in  C  n  :  h  ~  h  cos  4i, 

3 

in  £n  '  -2  sin  2i,  +  sin  4fc, 

4 

in  ns  6-8  cos  2l  +  2  cos  4  g. ; 


"t 

A4 

^4 

n4 


3  .  .  579  ,  „ 

64  +  4  cos  2t  -  -gj-  cos  4g, 

3  sin  2i  -  sin  4t, 

3  1  .  67 

^4  "  "3  cos  2i  +  cos  4fc, 


_3_ 

32 

_1_ 

32 

J_ 

64 

_5_ 

12 


i  -  Y2  sin  2)1  +  ff  sin  ^£» 

,  b4  *  0, 

1  .  67 

-  —  cos  2i  +  Ygj  cos  4i, 

29 

sin  2i  +  —  sin  4fc. 


Coefficient 

of 

in 

£2: 

in 

£n: 

2 

in 

n  : 

in 

c3= 

in 

£2n 

in 

£n2 

3 

in 

n  : 

4 

in 

£  : 

71  ,  37  31 

96  cos  £  ~  64  cos  3)1  “  192  cos  3)1 » 


37  ,  . 

48  Sin  1 


175  ,  .  41 

Y 92  sin  3£  +  192  sin  3£’ 


41  „  .  49  65 

-  cos  i  +  -32"  cos  3£  “  96  cos  3£’ 

15  „  .  27  » _  ,  3 

-  Jg  COS  i  +  Y2  cos  3£  +  32  COS  3£’ 


23 


101  .  13  ,  a  . 

-jg"  sin  i  ~  96  sin  3£  "  96  sin  3£’ 

91  „,73  37 

-  —  cos  i  +  64  cos  3i  -  YY2  cos  3£» 

65  .  ,  x  95  ,  31  , 

-  -J6  sin  £  +  37  sin  3£  ~  37  sin  3£» 


*S  cos  t  -  \  cos  3t  -  h  cos  5i, 
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in 

£3n: 

“  yy  sin  4  -  yy  sin  34  +  yy  sin  54, 

in 

cV: 

5  13  3 

~j~  cos  4  — y*  cos  34  +  y  cos  54, 

in 

Cn3: 

-  yy  sin  4  +  yy  sin  3  4  -  yy  sin  54, 

in 

4 

n  : 

2  cos  4-3  cos  34  +  cos  54, 

in 

£5: 

-  y  cos  4  -  yy-  cos  34  -  yy  cos  54, 

in 

C4n: 

13  1 

y  sin  4  +  y  sin  34  +  y  sin  54, 

in 

.3  2 
£  n  : 

-h  cos  4  +  h  cos  34  +  h.  cos  54, 

in 

£2n3: 

sin  4  +  *j  sin  34  -  ^  sin  54, 

in 

£n4: 

-2  cos  4+3  cos  34  -  cos  54, 

in 

n5: 

-20  sin  4-10  sin  34+2  sin  54; 

U5  ~  ^  cos  *■  +  “j^  cos  34  -  -^y^  cos  54, 

5  201  IRS 

Vj  =  -  yy  sin  4  +  -yy  sin  34  -  -yy-  sin  54, 

13  17 

£t-  =  -  yy  cos  34  +  yy  cos  54, 

n5  =  32  s*n  ^  -  <)(,  s^n  34  +  240  s^n 


A5 

-  0, 

B5 

0 

"  "  96 

9 

S 

5 

'  96 

cos 

4  -  11 
32 

COS 

34 

.  11 
+  48  ' 

n5 

5 

*  96 

sin 

4 

96 

sin 

34 

+  ^- 
240 

sin  54. 


